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Abstract 

We have made next to leading order QCD fit to the deep inelastic 
spin asymmetries on nucleons and we determined polarised quark and 
gluon densities. The functional form for such distributions was in- 
spired by the Martin, Roberts and Stirling fit for unpolarised case. 
In addition to usually used data points (averaged over x and Q 2 ) we 
have also considered the sample containing points with similar x and 
different Q 2 . It seems that splitting of quark densities into valence 
and sea contribution is strongly model dependent and only their sum 
(i.e., Au and Ad) can be precisely determined from the data. Inte- 
grated polarised gluon contribution, contrary to some expectations, is 
relatively small and the sign of it depends on the fact which sample 
of data points is used. 
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The final analysis of data on polarised deep inelastic scattering taken 
in E143 experiment at SLAC and SMC at CERN on proton and deuteron 
targets were published recently [|IJ, Together with the older data from 
SLAC @, BE, gi, CERN 0, |TT|, |T|, pi, g and DESY [15j one has 
quite a lot data on spin asymmetries. However, the newest data has smaller 
statistical errors and hence dominates in \ 2 fits. Study of polarised deep 
inelastic scattering were for the first time performed by SLAC- Yale group 
[0, f|. After the results of EMC group from CERN |IJ leading to the so 
called spin crisis there was enormous interest in studying polarised structure 
functions. It was suggested |16| that polarised gluons may be responsible for 
the little spin carried by quarks. Experimental groups have measured spin 
asymmetries on proton, neutron ( 3 He) and deuteron targets. On the other 
hand after calculation of two loop polarised splitting functions |17]] several 
next to leading (NLO) order QCD analysis were performed [|TJ, [T9, pT| 
making fits to the actually existing data and trying to determine polarised 
parton distributions. Determination of polarised gluon distribution was par- 
ticularly interesting in this context. The aim of this paper is similar. We 
want to perform next to leading order fit to the data taking into account the 
recently published data from the new analysis (SLAC E143 and CERN SMC 
experiments). We will divide the data into two groups. Many experimental 
groups published data |1|, |2j sets for the close values of x and different Q 2 
in addition to the "averaged" data where one averages over x and Q 2 (the 
errors are smaller and Q 2 dependence is smeared out). In principle when 
we take into account Q 2 evolution of polarised and unpolarised functions (in 
NLO analysis) the first group of data points, i.e. non averaged, should give 
better fit. In most of the fits to experimental data only second group of 
data namely with averaged x and Q 2 dependence was used. We will make 
fits to the both sets of data (the first group contains 374 points and the sec- 
ond 130 points) and compare them with the fits without Q 2 evolution taken 
into account (in other words assuming that asymmetries do not depend on 
Q 2 ). One should add that many experimental groups have not succeeded in 
finding Q 2 dependence for approximately the same value of x and different 
Q 2 m [13], [L4[]. The open question is if we really see the Q 2 evolution of 
structure functions from the existing experimental data. In our analysis we 
limit ourselves, as one usually does, to the data with Q 2 > 1 GeV 2 . As was 
already mentioned in our earlier papers [22| and was later stressed by other 
authors [18| making a fit to spin asymmetries and not directly to g\(x,Q 2 ) 
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we expect that the higher twist contributions are probably less important 
in such case. Experiments on unpolarised DIS provide information on the 
unpolarised quark densities q(x, Q 2 ) and G(x, Q 2 ) inside the nucleon. These 
densities can be expressed in term of q ± (x, Q 2 ) and G ± (x,Q 2 ), i.e. densities 
of quarks and gluons with helicity along or opposite to the helicity of the 
parent nucleon. The unpolarised quark densities are given by the sum of q + , 
q~ and G + , G~, namely: 

q = q + + q~, G = G + + G~ . (1) 

The polarised DIS experiments give also information about the so called 
polarised parton density, the difference of q + , q~ and G + , G~: 

Aq = q + -q~, AG = G + - G~ . (2) 

We will actually try to determine q ± (x, Q 2 ) and G ± (x, Q 2 ), in other words, 
we will have in some sense a simultaneous fit to unpolarised and polarised 
data. In principle the asymptotic x behaviour of q + and q~ will be taken from 
the unpolarised data. We will use fits of MRS (called R2) taking into 
account the behaviour at small x of quarks and gluon distributions obtained 
from experiments in Hera. It is of course very restrictive assumption that 
Aq and AG have the same behaviour (when the integrals of Aq and AG 
exist) as q and G. On the other hand the small x behaviour of unpolarised 
structure functions is determined from the x values of Hera much smaller 
than in the polarised case. We will also consider the integrals over the region 
measured in the experiments with polarised particles with the hope that in 
this case the behaviour of q ± and G^ could be more plausible. The values 
of integrals in the whole region (0 < x < 1) involving asymptotic behaviour 
taken from the unpolarised structure functions may be not as reliable as for 
the measured region. But it is an alternative to use of Regge type behaviour. 

It is known [I8| that the behaviour of the quark and gluon distributions 



in small x region is extremely important in extrapolation of calculation of 
integrals over whole < x < 1 range. It could happen that in Aq = q + — q~ 
(when we assume that q + , q~ and q = q + + q~ have similar x dependence) 
most singular x terms cancel (that is especially important in case of d va- 
lence quark, sea and gluon where the x behaviour is relatively singular). 
We will see how such description infers the fits and calculated parameters. 
With the less singular distributions for Ad v , AM (total sea polarisation) and 
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AG (gluon polarisation) there is no strong dependence of calculated quanti- 
ties in an unmeasured region but the fits have higher \ 2 ■ One of the main 
tasks of considering NLO evolution in Q 2 is the determination of the gluon 
contribution AG. In MS scheme AG(x, Q 2 ) comes in through the higher 
order corrections. In our fits we obtain AG relatively small contrary to some 
expectations. In addition when we use the averaged sample of data AG con- 
tribution is opposite to that in a set without averaging over x and Q 2 . In 
our fits the averaging over x and Q 2 (that have nothing to do with physics) 
changes the sign of the gluon contribution with rather small changes in quark 
region. It means that gluon contribution is extremely sensitive and can not 
be reliably determined. 

Let us start with the formulas for unpolarised quark parton distributions 
given (at Q 2 = 1 GeV 2 ) in one of the recent fits performed by Martin, Roberts 
and Stirling |2§. In this fit Ag 4 = 0.344 GeV and a s (M 2 z ) = 0.120. We 
have for the valence quarks: 

u v (x) = 2.251x~ a39 (l -x) 3 ' 54 (l -0.98^ + 6.51a;), 

d v (x) = 0.114a;- a76 (l-a;) 4 - 21 (l + 7.37^ + 29.9a;), (3) 

and for the antiquarks from the sea: 

2u(x) = 0.392M(V) -6{x), 

2d{x) = 0.392M(x) + 5(x), 

2s{x) = 0.196M(x), (4) 

2c(x) = 0.020M(x). 

In eq.(4) the singlet contribution M = 2[u + d + s + c] is: 

M{x) = 0.372T L15 (1 - s) 8 ' 27 (l + 1.13^/x + 14.4a;), (5) 

whereas the isovector 5 = d — u part: 

S(x) = 0.036aT L15 (l - a;) 8 - 27 (l + 64.9a;). (6) 

For the unpolarised gluon distribution we get: 

G{x) = 14.4x-°- 49 (l - x) 5 - 51 (l - 4.20^ + 6.47a;). (7) 



4 



We assume, in an analogy to the unpolarised case, that the polarised 
quark distributions are of the form: x a {l — x) 13 P 2 (\/x) , where P 2 (-Jx) is a 
second order polynomial in ^fx and the asymptotic behaviour for x — » and 
x — > 1 (i.e. the values of a and (5) are the same (except for AM, see a 
discussion below) as in unpolarised case. Our idea is to split the numerical 
constants (coefficients of P 2 polynomial) in eqs.(3, 5, 6 and 7) in two parts in 
such a manner that the distributions are positive defined. Our expressions 
for Aq{x) = q + (x) — q~(x) (q(x) = q + (x) + q~(x)) are: 

Au v (x) 
Ad v (x) 
AM(x) 

A5(x) 

AG(x) 

We will not consider AS, the parameter that breaks the isospin SU(2) 
symmetry of a sea, (we assume AS = 0) because one gets that the determi- 
nation of such parameters is not reliable. It is very important what assump- 
tions one makes about the sea contribution. From the MRS fit for unpolarised 
structure functions the natural assumption would be: As = Ad/2 = Au/2. 
If we add the condition that a$ value should be equal to the value deter- 
mined from the semileptonic hyperon decays, As is practically determined 
and though also nonstrange sea is determined. For comparison we will also 
consider such model. We assume that As (strange sea) is described by addi- 
tional parameters namely 

AM, = x- a65 (l - x) 8 - 27 (c u + c 2l Vx~) (9) 

where i = u,d,s and c Vl and c 2i for u and d are equal. For the strange quarks 
we have additional independent parameters. Comparing with the expression 
(5) we see that in AMj there is no term behaving like x~ 1A5 at small x (we 
assume that AMj and hence all sea distributions are integrable) which means 
that in AMj coefficient in front of this term have to be splitted into equal 
parts in AM 4 + and AM~ . We see that we have relatively strong singular 
behaviour of Ad v , AM and AG for small x values. For comparison we will 
also consider the model in which leading most singular terms are put equal to 
zero namely bi — Cu — d 1 — 0, that means that plus and minus components 



= x-°- 39 {l-xy- 5 \a 1 + a 2 Vx~ + a 3 x), 

= x-°- 76 (l-a;) 4 - 21 (6i + 62v^+63a;)> 

= x- a65 (l-x) 8 - 27 (c 1 + c 2V ^), (8) 

= x" a65 (l - x) 8 - 27 c 3 (l + 64.9a;), 

= x~ 0A9 (l-x) 5 - 51 c 3 (d 1 + d 2V ^+d 3 x). 



5 



have the same value for this powers of x that means that we investigate the 
dependence of the model on the leading x behaviour of Aq and AG (we know 
from the ref ]Tj3| that such dependence is strong). In the less singular models 



the dependence of calculated parameters in the unobserved region (below 
x < 0.003) is weak. In the earlier papers we considered the extrapolation of 
various calculated integrals below x = 0.003 up to assuming Regge type 
of behaviour for small x values. As will be discussed later the less singular 
models give slightly higher \ 2 . 

In order to get the unknown parameters in the expressions for polarised 
quark and gluon distributions we calculate the spin asymmetries starting 
from initial Q 2 = 1 GeV 2 for measured values of Q 2 and make a fit to the 
experimental data on spin asymmetries for proton, neutron and deuteron tar- 
gets. The asymmetry A±(x, Q 2 ) can be expressed via the polarised structure 
function Af(x,Q 2 ) as 

where R N = (FJ? — 2xF^)/2xF^ and F^ and F^ 1 are the unpolarised struc- 



ture functions. We will take the experimental from the ||24|| . Polarised 
structure function g±(x,Q 2 ) in the next to leading order QCD is related to 
the polarised quark, antiquark and gluon distributions Aq(x, Q 2 ), Aq(x, Q 2 ), 
AG(x, Q 2 ), in the following way: 

2 - 2tt 

+ Aq(x,Q 2 )) + j5c g *AG(x,Q 2 )]} (11) 

with the convolution * defined by: 

(C*q)( Xl Q 2 )= f -C(-)q( Zl Q 2 ) (12) 

Jx Z Z 

The explicit form of the appropriate spin dependent Wilson coefficient in 



the MS scheme can be found for example in ref. []I7[] • The NLO expressions 
for the unpolarised (spin averaged) structure function F N (x,Q 2 ) is similar 
to the one in eq.(ll) with Aq(x, Q 2 ) — > q(x, Q 2 ) and the unpolarised Wilson 
coefficients are given for example in BR EBfl. 
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The Q 2 evolution of the parton densities is governed by the DGLAP 
equations [pTfl . For calculating the NLO evolution of the spin dependent 
parton distributions Aq(x, Q 2 ) and AG(x, Q 2 ) and spin averaged q(x, Q 2 ) 
and G(x, Q 2 ) we will follow the method described in p6| [L8[]. We will cal- 
culate Mellin n-th moment of parton distributions Aq(x, Q 2 ) and AG(x, Q 2 ) 
according to 

Aq n {Q 2 ) = [ dxx^Aqix.Q 2 ) (13) 
Jo 

and then use NLO solutions in Mellin n-moment space in order to calculate 
evolution in Q 2 for non-singlet and singlet i.e. of AS n (Q 2 ) = J2 q [Aq n (Q 2 ) + 
Aq n {Q 2 )} and AG(Q 2 ). 

In calculating evolution of S(Q 2 ) and AG(Q 2 ) with Q 2 we have mixing 
governed by the anomalous dimension 2x2 matrix. We used explicit formulae 
given in ||26|| . Having evolved moments one can insert them into the n-th 
moment of eq.(ll). 

9i(Q 2 ) = \lle 2 q {Aq n {Q 2 ) + Aq\Q 2 ) + ^{5c n q -{Aq\Q 2 ) 

+ A^(Q 2 )) + j5c n g -AG n (Q 2 )]} (14) 

and then numerically Mellin invert the whole expression. In this way we get 
gi (x,Q 2 ). The same procedure is applied taking the appropriate formulas 
giving the different Q 2 dependence and the correction coefficients for the un- 
polarised structure functions. Having calculated the asymmetries according 
to equation (10) for the measured in experiments value of Q 2 we can make a 
fit to a measured asymmetries on proton neutron and deuteron targets. We 
will take into account for proton 7 points of E80 0] and 16 points of E130 [§J 



of SLAC experiments, 10 points of EMC [l(J and 59 points of SMC § from 
CERN and 82 points of E143 |]J from SLAC. For deuteron we have 65 points 
from SMC and 82 points from E143 |1] whereas for neutron 33 points of 
E142 |J and 11 points of E154 M experiments from SLAC and 9 points from 



DESY Hermes experiment []To| . The last two sets of data from E154 and Her- 
mes are taken in order to have more data from neutron target and to balance 
huge number of points from proton and deuteron targets. All together we 
have 374 data points and together with the assumed g$ = 0.579 ±0.1 value 
we have 375 points (174 for proton, 147 for deuteron and 53 for neutron). 
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<2i = 


0.66, 


a 2 = 


-4.21, 


«3 


= 14.6, 


h = 


-0.02, 


b 2 = 


-0.84, 


h 


= -1.74, 


Clu — 


c ld = -0.28, 


C2u — 


c 2d = 3.08, 






Cls = 


-0.42, 


C2s — 


-1.15, 


c 3 


= (input), 


di = 


2.201 


d 2 = 


-22.47 


dz 


= 42.20. 



We get the following values of parameters from the fit to all existing 
(above mentioned) data for Q 2 > lGeV 2 for spin asymmetries: 



(15) 



The resulting y 2 per degree of freedom is y 2 /N^f = Jfjpi3~ = 0.853. 

The obtained quark and gluon distributions lead for Q 2 =1 GeV 2 to 
the following integrated quantities: Am = 0.77 (Au v = 0.70, 2Au = 0.07), 
Ad = -0.49 (Ad v = -0.56, 2Ad = 0.07), As = -0.15. These numbers give 
the following predictions: AS = 0.13 , AM = 0.00, AG = 0.22, T\ = 0.113, 

= -0.062, Tf = 0.024, g A = Au - Ad = 1.26. 

We have relatively small positively polarised sea for up and down quarks 
and stronger negatively polarised sea for strange quarks. The gluon polar- 
isation is positive but very small. The value of g A was not assumed as an 
input in the fit and comes out correctly. 



As was already stressed in [[31|, the asymptotic behaviour at small x of 
our polarised quark distributions is determined by the unpolarised ones and 
these do not have the expected theoretically Regge type behaviour or pQCD 
which is also used by experimental groups, to extrapolate results to small 
values of x. Some of the quantities change rapidly for x < 0.003. 

Now we will present quantities integrated over the region from £=0.003 
to x=l (it is practically integration over the region which is covered by the 
experimental data, except of non controversial extrapolation for highest x). 
The corresponding quantities are Am = 0.78 (Au v = 0.67, 2An = 0.11), 
Ad = -0.42 (Ad v = -0.53, 2AJ = 0.11), As = -0.12, AS = 0.24, AM = 
0.11, AG = 0.06. We can also calculate T p , T n and T d in the measured 
region for Q 2 = 5 GeV 2 and compare them with the quantities given by the 
experimental groups. 

We get rf = 0.119, = -0.078 and Tf = 0.019 in the whole region, 
whereas in the region between x = 0.003 and x = 0.8 (covered by the data) we 
have rf = 0.125, r™ = -0.051 and Tf = 0.034. The experimental group SMC 



present |20| the following values in the measured region (for Q 2 = 5 GeV 2 ): 

T{ = 0.130 ±0.007, 
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1 1 
1 1 



-0.054 ±0.009, 
0.036 ±0.005. 



(16) 



The world average for such Q 2 is [20] for the whole region (0 < x < 1): 

r? = 0.121 ±0.018, 

= -0.075 ±0.021, (17) 
r? = 0.021 ±0.017. 

Our results are in good agreement with given experimental values. For com- 
parison we have also made fits using formulas of the simple parton model 
(as in our papers before f22|) neglecting evolution of parton densities with 
Q 2 . More detailed result of these fits (integrated densities and so on) will be 
given later. 





Figure 1: Plots of proton spin asymmetry predicted by our basic fit (for experimental Q 2 ). 
The data points from different experiments (E14-3, SMC) with total errors are also shown. 

In figs.l, 2 and 3 we present the comparison of our basic fit with measured 
asymmetries for proton (1), deuteron (2) and neutron (3) targets. The curves 
are obtained by joining the calculated values of asymmetries corresponding 
to actual values of x and Q 2 for measured data points. For the same value 
of x we have experimental points corresponding to different Q 2 values. We 
see that distribution of experimental points is much bigger than the lengths 
of vertical lines measuring the changes of influence of evolution in Q 2 for 
different values of Q 2 and the same value of x. It seems that with such errors 
it is difficult to see Q 2 dependence of asymmetries. 
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Figure 2: Plots of deuteron spin asymmetry predicted by our basic fit (for experimen- 
tal Q 2 ). The data points from different experiments (ElJj.3, SMC) with total errors are 
also shown. 
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Figure 3: Plots of neutron spin asymmetry predicted by our basic fit (for experimen- 
tal Q 2 ). The data points from different experiments (E14-2, E154) with total errors are 
also shown. 

For asymmetries the curves with Q 2 evolution taken into account and 
evolution completely neglected do not differ very much so we do not present 
them. We see that in the case of g± function for proton (fig.4) the dashed line 
corresponding to the fit with no evolution in Q 2 taken into account (parton 
model) that for small x is continuation of E143 data (small errors) but lays 
within experimental errors of SMC results. May be it could be considered as 
some tendency in evidence for seeing Q 2 dependence in the data but certainly 
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not a very strong one. In the deuteron and neutron data the effect is even 
less pronounced. 




X 



Figure 4: Plots of structure function gi(x, Q 2 ) obtained in our fit (for experimental values 
of Q 2 ) and compared with data points for proton, deuteron and neutron targets. The solid 
curves are predictions for SMC or E142 experimental points, whereas dot-dashed ones are 
for El 43 or El 54 data. The dashed curves are predictions from the fit with no Q 2 evolution 
of considered structure functions. 

Polarised quark distributions for up and down valence quarks as well as 
non strange, strange quarks and gluons for Q 2 = 1 GeV 2 are presented in 
figure 5. They are compared with the distributions when Q 2 evolution is not 
taken into account and with the corresponding unpolarised quark and gluon 
distributions. We see that especially in the case of polarised gluon distri- 
bution function does not resemble the distribution function for unpolarised 
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case. This function is also quite different from the gluon distribution (given 
in |8f) used to estimate AG/G in COMPASS experiment planned at CERN 

im- 

Fixing the value of g 8 is very important for the fit. When we relax the 
condition for g 8 = 0.579, x 2 g° es a little bit down to the value 308.65. We 
get the fit with the parameters not very different from our basic fit but with 
very small g 8 (0.03) and bigger AS = 0.27, (AS = 0.37 for 0.003 < x < 1) 
and positive As = 0.08, Am = Ad = 0.03. It means that fixing the value of 
g 8 equal to experimental value (gotten from hyperon /3-decays data) enforces 
the negative value of As between -0.15 and -0.12. The obtained solution 
without fixing g 8 value is somehow dual to our fit, g 8 very small and AS 
relatively large comparing with g 8 close to its experimental value and AS 
rather small. 

In order to check how the fit depends on the assumptions made about the 
sea contribution we have also made fit with Am = Ad = 2As, the assumption 
that follows directly from MRS unpolarised fit. The \ 2 value per degree 
of freedom \ 2 I^df = J^rfr = 0.856 is a little bit worse. In this case we 
have Au = 0.76 (Am, = 0.98, 2Au = -0.21), Ad = -0.48 (Ad v = -0.26, 
2Ad = -0.21), As = -0.11, AS = 0.18, AM = -0.53, AG = 0.28. Hence, 
we see that with the different assumption about sea behaviour the overall sea 
contribution changes quite drastically. The quantity As must be negative 
in order to get experimental value for g 8 and because of our assumption 
(Am = Ad = 2As) we obtain relatively big negative values of non strange 
sea for up and down quarks. We see that the values for sea polarisation 
depend very strongly on the taken assumptions (in many papers Jl^, |19|, p8| 
SU(3) symmetric sea is assumed that also together with fixing of g 8 value 
gives negative non strange sea). It means that the sea contribution is not very 
well determined. On the other hand Am = Am„ + 2Au and Ad = Ad v + 2 Ad 
practically do not change (however, Am„ and Ad v also change). Also AG 
does not change. 

Looking at the dependence of unpolarised quark and gluon densities we 
see that the most singular behaviour for small x we have for d v (x), M(x) 
and G(x). For comparison we have investigated the model when in polarised 
densities these most singular contributions are absent. In this case Ad v , 
AM and AG are ^/x less singular than in our basic fit. For such a fit we get 
X 2 /Ndf = 375^9 = 0.864 i.e only slightly higher than in our basic fit. We get in 
this case: Am = 0.78 (Au v = 0.71, 2Am = 0.07), Ad = -0.41 (Ad v = -0.48, 
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Figure 5: Polarised quark (xAu v , xAd v ), antiquark (xAu, xAs) and gluon distributions 
(xAG) predicted by our basic jit at Q 2 = IGeV 2 (solid curve). For comparison prediction 
for such quantities for the fit without Q 2 evolution taken into account (long- dashed curves) 
as well as nonpolarised quark, antiquark and gluon distributions from J13| / (dashed ones) 
are also shown. 
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2Ad = 0.07), As = -0.10, AS = 0.27, AM = 0.05, AG = -0.40. If we do 
not modify AG(x, Q 2 ) omitting the most singular term AG remains positive. 
In this fit integrated quantities taken over the whole range of < x < 1 and 
in the truncated one (0.003 < x < 1) differ very little. AG is negative 
opposite to the basic fit. So it is possible to get the fit of comparable quality 
to our basic fit with practically no change of integrated quantities in the 
region between x = and x = 0.003. For Q 2 =l GeV 2 we have = 0.121 
and = -0.044 (to be compared with = -0.062) and relatively big 
AS = 0.27. 

The obtained results can be compared with the fit when instead of 374 
points for different x and Q 2 values we take spin asymmetries for only 130 
data points with the averaged values for the same x, averaged Q 2 and smaller 
errors. We have then for proton target points obtained at CERN (10 from 
EMC, 12 for SMC) and at SLAC (4 from E80, 8 from E130 and 28 for E143). 
For deuteron we take into account 12 points from SMC and 28 from E143 
whereas for neutron target data points from SLAC (8 from E142 and 11 from 
E154) and DESY(9 from Hermes). In this fit first of all the errors are smaller 
than in our basic fit and the ratio of number of neutron to number of deuteron 
and proton data points are different. It seems that the influence of neutron 
points is stronger than in basic fit x 2 /Ndf — if?Zi3 = 0.844 is a little bit better 
than in our basic fit. The integrated values for quark and gluon densities 
are:_ Am = 0.80 (Au v = 0.57, 2Au = 0.23), Ad = -0.54 (Ad v = -0.78, 
2Ad = 0.23), As = -0.16, AS = 0.10, AM = 0.31, AG = -0.69 and 
qa = 1.35. We see that averaging over x and Q 2 and different numbers of 
data points leads to not a very different fit (for example \ 2 = 101.2 when 
the parameters of the basic fit are used) but the values for integrated valence 
densities and nonstrage sea contribution are different (Au and Ad do not 
differ very much and the difference is smaller for integrated quantities in the 
region 0.003 < x < 1). Integrated gluon density is negative. The fits become 
more similar when in the average fit we fix the qa value, i.e. the condition 
for valence contribution. 

The fact that the averaged and non averaged samples of data points 
results for valence quark densities and sea contributions are different and in 
the case of integrated gluon density even sign is different means that these 
quantities are not very precisely determined in the fits. After assuming g 8 
value from experiment and fixing As contribution we can determine Au and 
Ad values (the splitting of Aq in Aq v and 2Ag depends on data sample and 
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assumptions about sea contribution). It could be that differences that come 
out in comparing fits to the averaged and non averaged data are connected 
with the fact that rather singular polarised parton distributions are able to 
pick up differences in two experimental data samples, due to different number 
of neutron to proton and deuteron data points. The fit with less singular 
densities (Ad v , Am, Ad and AG) with x 2 a little bit higher 104.87 {x 2 /N DF 
= T§T^9 ) =0.860) is for averaged data points nearly identical to that for non 
averaged one. AG is like in previous case small and negative AG = —0.4. 

As was already mentioned before we have also made for comparison fits 
neglecting evolution of parton densities with Q 2 (formulas from the simple 
parton model). We get for non averaged data sample x 2 I^df = J^pfj)) = 
0.872 (higher than in our basic fit where: x 2 /^df = 0.853): Am = 0.68 
(Au v = 0.46, 2Au = 0.22), Ad = -0.41 (A^ = -0.63, 2Ad = 0.22), 
As = -0.16, AS = 0.11, AM = 0.27, r? = 0.120, = -0.062. For 
averaged data points we get x 2 I^df = 1 1 3 1 3 1 7 L 9 ) = 0.858 (this number should 
be compared with x 2 I^df = 0.844, the corresponding quantity from NLO 
fit) and we have: Au = 0.69 (Au v = 0.40, 2Au = 0.29), Ad = -0.42 
(Ad v = -0.72, 2AJ = 0.29), As = -0.16, AS = 0.11, AM = 0.42. Hence, 
X 2 is smaller in the case of averaged sample. We see that also in this case 
the value Am and Ad and As are practically the same and there are some 
shifts in valence values and non strange sea contribution (the differences are 
smaller in the 0.003 < x < 1 region). Similarly to the case with evolution 
taken into account the fits are nearly identical when we make them with less 
singular Ad v , Au = Ad and AG densities (x 2 is in this case higher than in 
the basic fit for non averaged data sample or for averaged one). 

It has been pointed out that the positivity conditions could be re- 
strictive and influence the contribution of polarised gluons. We have also 
made a fit to experimental data without such assumption for polarised par- 
tons. The x 2 value does not changed much x 2 I^df = Jffzf| = 0.852 and 
we get Am = 0.79 {Au v = 0.35, 2Am = 0.44), Ad = -0.50 (Ad v = -0.94, 
2AJ = 0.44), As = -0.12, AS = 0.17, AM = 0.76, AG = 0.12. Hence, 
sea contribution for non strange quarks became big and there are shifts be- 
tween valence and sea contribution. The gluon contribution does not change 
staying close to zero for non averaged data and small and negative for av- 
eraged data. The most important restriction is for the valence d quark (in 
the case without positivity assumption it becomes big close to -1). Relax- 
ing the positivity of other quantities does not change the fit. It seems that 
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without positivity the values of parton distributions are strongly pushed in 
the direction of small contribution of valence u quark and big (and negative) 
valence d quark. The gluon contribution is not modified very much. Hence, 
we decided to use the positivity assumptions in our fit. 

In summary we have made fits to two samples of data with averaged 
x and Q 2 values and nonaverged ones (adding averaged neutron data from 
E154 and Hermes experiments). To check the influence of different model 
dependent assumptions we consider fits without fixing g$ value, with mod- 
ified sea contribution and less singular behaviour for valence d quark, sea 
contribution and gluon densities. For comparison we have also considered 
fits to the simple parton model neglecting Q 2 dependence of parton densi- 
ties. It seems that splitting of integrated densities Au, Ad in valence and sea 
contribution is model dependent (Au and Ad were more or less the same). 
The integrated gluon contribution comes out relatively small and of different 
sign in averaged and nonaveraged data sample. The comparison of g\ for 
fits without evolution and with Q 2 evolution of parton densities taken into 
account is given. It seems that experimental accuracy is still not enough to 
make precise statements about polarised quark and gluon densities. 
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